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We study the role the proton-electron gas interaction has on the formation of nuclear structures 
in neutron star crusts. Using a classical molecular dynamics model we study isospin symmetric and 
asymmetric matter at subsaturation densities and low temperatures varying the Coulomb interaction 
strength. The effect of such variation is quantified on the fragment size multiplicity, the inter-particle 
distance, the isospin content of the clusters, the nucleon mobility and cluster persistence, and on 
the nuclear structure shapes. We find that the proton-electron Coulomb interaction distributes 
matter more evenly, disrupts the formation of larger objects, reduces the isospin content, modifies 
the nucleon average displacement, but does not affect the inter-nucleon distance in clusters . The 
nuclear structures are also found to change shapes by different degrees depending on their isospin 
content, temperature and density. 
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I. INTRODUCTION 

The raison de vivre of neutron stars is a balance be- 
tween nuclear, Coulomb and gravitational forces. In vac- 
uum, the weak decay n —> p+e~ + v e is exothermic with a 
net energy release of 0.778 MeV as kinetic energies of the 
emitted particles. In conditions of thermal and chemical 
equilibrium, the opposite process, p + e — > n + v, which 
is endothermic, would work to restore the original ratio 
of protons and neutrons, this, however, is not the case at 
neutron star pressures and densities. Due to Pauli prin- 
ciple arguments, the neutron decay is inhibited and such 
processes yield to an accumulation of neutrons; the neu- 
tron dominance is expected to start in the stellar crust 
when Z/A ~ 0.13 jlfl, while the outer layers maintain 
larger neutron to proton ratios in interesting structures 
dictated by the varying percentage of protons and neu- 
trons in an all-embedding electron gas. 

The electron gas is expected to play an important role 
in the formation of the nuclear structures. As the den- 
sity increases as a function of depth in the crust, the 
electrons will transform protons into neutrons by inverse 
j3 decay making the neutron-to-proton ratio more asym- 
metric and forming neutron-rich nuclei. In addition to 
being responsible for producing neutrons, the electron 
cloud also creates a Coulomb potential that surrounds 
the rest of the matter contributing to define its shape. 

At densities and temperatures expected to exist in 
neutron star crusts (p < po and T < 1.0 MeV, with 
Po denoting the normal nuclear density), nucleons form 
structures that are substantially different from the quasi- 
spherical nuclei we are familiar with in our part of the 
universe. Such structures, which have been dubbed "nu- 
clear pas ta" , have been investigated using various mod- 
els [3-[l3l] which have shown them to be the result of the 
interplay of nuclear and Coulomb forces in an infinite 



medium. The structure of the nuclear pasta is expected 
to play an important role in the study of neutrino opac- 
ity in neutron stars [idj . neutron star quakes and pulsar 
glitches [l|. 

In a previous study, a combination of classical molec- 
ular dynamics and fragment recognition algorithms (l3j 
was shown to be very effective in the study of the pasta 
structures. By using a microscopic view of the nucle- 
ons, it was found that a combination of topological tools 
can be used to classify the different structures into rec- 
ognizable patterns; this opens the door to cross-model 
comparison between structures obtained with different 
approaches and to a quantifiable analysis of the effect 
the nuclear and Coulomb forces have on the pasta for- 
mation. This investigation will focus on the role of the 
Coulomb forces on the creation of the pasta. 

The presence of the electron gas is bound to have, at 
least, a two-fold effect. On one side, its interaction with 
the protons should certainly modify their positioning in 
space as well as their kinetic energy and, on another 
hand, the fact that the gas fills all openings should re- 
duce the long range proton-proton interactions with an 
screening-like effect . 

In a 2003 investigation [l2|, the screening effect of an 
electron gas on cold nuclear structures was investigated 
using a static liquid-drop model, and it was found that 
main effect of the gas screening was to extend the range of 
densities where bubbles and clusters appear. Indeed, the 
screening was found to be of minor importance although 
the study, being static, did not include any spatial and 
dynamical effects the gas should have produced on the 
nucleons. 

A 2005 study pj] used a density functional method 
to investigate charge screening on nuclear structures at 
subnuclcar densities but still at zero temperature; in par- 
ticular, cases with and without screening were directly 
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compared. The main results of the study were to ex- 
tract the nucleon density profiles and to quantify the 
spatial rearrangement of the proton and electron charge 
densities. Once again, it was found that the density re- 
gion in which the pasta exists becomes broader when the 
Coulomb screening is taken into account, mainly due to 
the rearrangement of the protons; the authors remark 
the importance of extending such study to finite temper- 
atures and with dynamical models. 

Recently, in an extension of previous works flil Il7l. [l8j . 
Pickarewicz and Toledo Sanchez performed Monte 
Carlo simulations directly on 5,000 nucleons to deter- 
mine the lowest proton fractions that arc compatible 
with /3 equilibrium in neutron star crust environments; 
a feature that depends finely on the Coulomb interac- 
tion between the protons and electron gas. The study, 
which approximated the Coulomb interaction via an 
Ewald summation, was performed at a single tempera- 
ture (T = 1.0 MeV) and it succeeded in finding a rela- 
tionship between the proton fractions and the density of 
the electronic contribution. In closing, the authors rec- 
ommended that further investigations be performed with 
molecular dynamics simulations to study dynamical ob- 
scrvables and to examine the role of the temperature in 
the study of stellar crusts. 

Thus, the motivation of this work is to study the effect 
the proton-electron gas interaction has on the formation 
of the nuclear pasta using a dynamical model, at non-zero 
temperatures, and in a systematic way that takes advan- 
tage of the tools developed in our previous work fl3| . 
In the next section the model used (classical molecular 
dynamics) is briefly introduced along with the cluster 
recognition algorithms and analysis tools employed. Sec- 
tion HIT] presents a series of results of the simulations and 
analyses performed. We summarize some of our find- 
ings in Section IIVI and present a quick look at work in 
progress. 



II. MODEL USED 

The model used here was developed to study nuclear 
reactions from a microscopic point of view, and it is 
composed of a molecular dynamics code retrofitted with 
cluster recognition algorithms and a plethora of analy- 
sis tools. The justification for using this model in stellar 
crust environments was presented elsewhere fl3j , here we 
simply mention some basic ingredients of the model. 

The classical molecular dynamics model (CMD), first 
introduced in (20| . has been used in heavy- ion reaction 
studies to explain experimental data [2ll|, identify phase 
transitions and other critical phenomena I22l425l | , explore 
the caloric curve of nuclear matter (2(| |27| and isoscal- 
ing [2^, . Synoptically, CMD uses two two-body po- 
tentials to describe the motion of "nucleons" by solving 
their classical equations of motion. The potentials were 



FIG. 1: (Color online) Characteristic structures obtained with 
(a=l) and without (a=0) the proton-electron gas interaction 
for symmetric matter (x = 0.5) at density p = 0.015 fm~ 3 
and temperature T = 0.1 MeV. 

developed phenomenologically by Pandharipande [30^ : 

V np (r) = V r [exp(-iJ, r r)/r - exp(-n r r c )/r c ] 

- V a [exp(-fi a r)/r - exp(-p, a r c )/r c ] 
V NN (r) = V [exp{-fior)/r -exp(-fior c )/r c ] , 

where V np is the potential between a neutron and a pro- 
ton, and Vnn is the repulsive interaction between either 
nn or pp. The cutoff radius is r c = 5.4 fm and for r > r c 
both potentials are set to zero. The Yukawa parame- 
ters n r , [i a and /zo were determined by Pandharipande to 
yield an equilibrium density of po = 0.16 /m -3 , a binding 
energy E(p ) = 16 MeV/nucleon and a selectable com- 
pressibility of 250 MeV for the "Medium" model, and 
535 MeV for the "Stiff"' 0. We remark that, since 
the potential Vnn does not permit bound states between 
identical nucleons, pure neutron matter is unbound and, 
at a difference from potentials used by other models [3] , 
the Pandharipande potentials have a hard core. 

The main advantage of the CMD model is the possi- 
bility of knowing the position and momentum of all par- 
ticles at all times as this allows the study of the structure 
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of the nuclear medium from a microscopic point of view. 
The outcome of CMD, namely, the time evolution of 
the particles in (r, p), can be used as input in anyone of 
the several cluster recognition algorithms that some of us 
have designed for the study of nuclear reactions [3lTj33l | . 

As explained elsewhere [13l l3~il |35| , the lack of quan- 
tum effects, such as Pauli blocking, -perhaps the only 
serious caveat in classical models- ceases to be relevant 
in conditions of high density and temperature (such as 
in heavy-ion reactions) or in the low-density and low- 
temperature stellar environments. 

To simulate an infinite medium, systems with thou- 
sands of nuclcons were constructed using CMD under pe- 
riodic boundary conditions. Cases symmetric in isospin 
(i.e. with x = Z/A = 0.5, 1000 protons and 1000 neu- 
trons), and cases with asymmetric isospin (with x = 0.3 
comprised of 1000 protons and 2000 neutrons) were con- 
structed in cubical boxes with sizes adjusted to have den- 
sities between p = 0.01 fm < p < po; our previous 
study, flo| , presented an smorgdsboard of structures ob- 
tained through this method. 



A. Coulomb interaction 

The neutron star crust is embedded in a degener- 
ate electron gas. Although the nucleon-electron sys- 
tem is overall neutral and /3-equilibrated, to take into 
account the Coulomb interaction (which is of infinite 
range) it is necessary to use some approximation. The 
two most common approaches are the Thomas-Fermi 
screened Coulomb potential (used, e.g., in QMD and 
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FIG. 2: Effect of the screened potential on the fragment size 
multiplicity. Plots show the distribution of cluster sizes ob- 
served in 200 configurations of x — 0.3 nuclear matter at den- 
sity p = 0.015 fm' 3 and temperatures T = 1.0 MeV (top) 
and T = 0.1 MeV (bottom). The figure on the left panel 
correspond to configurations with Coulomb interaction, and 
those on the right to the case without such potential. 
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FIG. 3: Examples of the radial correlation function for vary- 
ing strengths of the Coulomb potential: a = 1 (full Coulomb), 
0.8, 0.2, and (without Coulomb). 



the Ewald summation procedure 1361; CMD has been 
used with both approximation (see 1 1 31 ] for a comparison 
of methods under CMD) but for computational reasons 
we opt for the use of a screened Coulomb potential. 

We approximate the electron gas as a uniform ideal 
Fermi gas with the same number density as the pro- 
tons. Its inclusion in the CMD is through the use of 
a screened Coulomb potential obtained from the Pois- 
son equation, V^ Scr \r) = (e 2 /r) exp(— r/A), where the 
screening length, A, is related to the electron mass, Fermi 
momentum, and number density. Here, pragmatically, 
and after testing various lengths, we use the prescription 
used before [1,(3 an d se t A = 10 fm, which satisfies the 
requirement of being sufficiently smaller than the size of 

1/3 

the simulation cell, L = (A/p) 1 . 



B. Simulation procedure 

The trajectories of the nucleons arc then governed 
by the Pandharipande and the screened Coulomb po- 
tentials. The equations of motion are solved using a 
symplectic Verlet integration with energy conservation of 
0(0.01%). The nuclear system is force-heated or cooled 
using isothermal molecular dynamics with the Ander- 
sen thermostat procedure [37| with small temperature 
steps. To simulate stellar crust conditions we operate 
in the range of 0.1 < T < 1.0 MeV and densities 
p = 0.01 fm < p < po; in comparison to heavy-ion 
reactions, nucleons in neutron star crusts are practically 
in a semi-frozen state. 

To obtain reliable statistics we sample each T, p and x 
configuration 200 times in an stationary ergodic process. 
That is, after reaching equilibrium at a given set of values 
of T, p and x, the (r, p) values of all nucleons are stored 
for further processing while the system continues its evo- 
lution for a "long" time (4000 time steps) until it reaches 



4 




ol , , i , , i , , , I o 

20 40 60 80 100 120 140 160 180 200 

Configuration Number 



0.1 1 1 1 1 1 

250 500 750 1000 

Fragment Size 



FIG. 4: x content of the clusters formed in 200 configurations 
of asymmetric matter (x = 0.3) at density p = 0.015 fm~ 3 
and temperature T — 1.0 MeV. The top panel shows the 
case with Coulomb and the bottom one the case without such 
interaction. Notice that the abscissas have different scales. 



FIG. 5: RMS displacement of nucleons during the evolution 
of the stationary ergodic process with Coulomb (continuous 
line) and without Coulomb (dashed line) as labeled; see text 
for details. Notice that at low temperatures the mobility with 
Coulomb is smaller than without Coulomb, while at higher 
temperatures the effect is the opposite. 



C. Analysis tools 



a state independent of the previous one, point at which 
the procedure is repeated for a total of 200 times. The 
recorded nucleon positions and momenta are then used 
later to identify clusters and to characterize the structure 
as explained next. 

Given that, for the densities and temperatures of in- 
terest, the nucleons do not have large mobilities, the 
task of identifying self-bound cumuli is sufficiently sim- 
ple as to be achieved with a basic "Minimum Spanning 
Tree" (A/ ST) algorithm (HEIl- The MST finds clus- 
ter membership simply by identifying those particles that 
are closer to each other than some clusterization radius, 
which here is set to 3.0 fm based on the range of the 
nuclear potential and densities investigated; the generic 
method was modified to recognize fragments that extend 
into any of the adjacent periodic cells. 



Once the fragments have been identified, their global 
properties are studied as follows. MST yields infor- 
mation about the fragment multiplicity; in Ref. fl3j 
systems at x = 0.3 and 0.5 and at various densities 
(0.01 /m~ 3 < p < 0.015 fm~ 3 ) and temperatures 
(0.1 MeV < T < 1.0 MeV) exhibited increases in mul- 
tiplicity with temperature (in opposition to the rise-fall 
that takes place in intermediate-energy heavy-ion reac- 
tions) as well as the formation of percolating structures 
that encompass most of the nucleons in the simulation. 

At a microscopic level, the dynamics of the nucleons 
in such systems can be quantified through their average 
displacement as a function of "time" , i.e. through the 
time steps of the simulation; for the cases mentioned be- 
fore (cf. [l3[) the mobility between configurations was 
found to be between 2 and 20 fm. Likewise, the mi- 
croscopic stability of the clusters can be gauged through 
the "persistency" [HI, [H[ which measures the tendency 
of members of a given cluster to remain in the same 
cluster; for the cases previously studied, colder clusters 
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(T < 0.5 MeV) exchange less than half of their nuclcons 
during 200 ergodic configurations, while those at higher 
temperatures (T > 0.8 MeV) replace up to 90% of their 
constituents. 

Another interesting descriptor is the isospin content of 
each cluster produced, i.e. their x values; in the pre- 
vious study it was found that for the case of x = 0.3, 
small clusters (A < 10) tend to have less protons than 
the global average, while such effect was not present in 
the case of x = 0.5. Other global characterization tool is 
the pair correlation function, g(r), which gives informa- 
tion about the spatial ordering of the nuclear medium by 
comparing the average local density to the global den- 
sity; in the previous study g(r) showed that nuclcons in 
clusters have an interparticlc distance of about 1.7 fm at 
all studied densities, and was useful in determining the 
onset of lasagna-like structures. 





x=0.3, T=1.0 MeV, p=0.015 fm" 3 
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FIG. 6: Persistence with and without Coulomb as it evolves 
through 200 configurations of the stationary ergodic process. 

Proceeding beyond global measures, the shapes of nu- 
clear structures can be characterized by their volume, 
surface area, mean curvature, plus an interesting con- 
struct known as the Euler characteristic; these four ob- 
jects comprise what is known as the "Minkowsky func- 
tionals" which completely describe all topological prop- 
erties of any three-dimensional object. The evaluation 
of the mean curvature (averaged over the whole surface 
of the object) and Euler number (a topological invariant 
connected to the vertices, edges and faces of polyhedra), 
however, requires the mapping of the nuclear clusters 
into a polyhedra; procedure that can be accomplished 
through the Michielsen-De Raed algorithm [ll], [39| . 

In [l3j it was shown that generic structures, such 
as "gnocchi", "spaghetti", "lasagna" and "crossed- 
lasagnas" or "jungle gym" and their inverse structures 
(with voids replacing particles and viceversa), all have 
well defined and distinct values of the mean curvature 
and Euler numbers with magnitudes dictated by the over- 
all size of the structure, i.e. by the number of particles 
used. It should be remarked that the assignment of cur- 
vature and Euler values to a given structure becomes 



relative to the overall scale. This becomes important 
for structures with near-zero Euler number which sig- 
nal spaghetti, lasagna and their anti-structures; curled 
spaghetti rods of a given length, for instance, can easily 
oscillate between positive and negative near zero Euler 
numbers by curling or uncurling a segment. This classifi- 
cation, which was respected at all x values, densities and 
temperatures studied, is shown in Table HI 

It is with this model and arsenal of tools that we now 
embark on a study of the effect the proton-electron gas 
interaction has on the formation of the nuclear pasta. 



TABLE I: Classification Curvature - Euler 





Curvature < 


Curvature ^ 





Curvature > 


Euler > 
Euler ~ 
Euler < 


Anti- Gnocchi 
Anti- Spaghetti 
Anti- Jungle Gym 


Lasagna 


Gnocchi 
Spaghetti 
Jungle Gym 



III. THE EFFECT OF COULOMB 

Along the line of previous studies [HI, [l|| [lj| and ad- 
hering to their advise, in this study we use a molecular 
dynamics model to study spatial and dynamical effects 
the electron gas may produce on the nuclear structures 
at non-zero temperatures. Specifically we extend the ap- 
proach of one of the previous studies and compare pasta 
structures obtained with and without the proton-electron 
gas interaction as well as with varying strengths of it. 

In summary, pasta structures were obtained at sub- 
saturation densities (0.015 fm~ 3 < p < 0.072 /to" 3 ) 
and low temperatures (0.1 MeV < T < 1.0 MeV) 
for symmetric (x = 0.5) and asymmetric (x = 0.3) 
cases; altogether 200 simulations were carried out per 
each combination of {p, T, x}. To see the effect of the 
Coulomb interaction, each case was "cooked" repeatedly 
with a screened Coulomb potential with a varying am- 
plitude, i.e. with aV^ Scr \ where < a < 1. For a fair 
comparison, all corresponding cases with different val- 
ues of a were produced with identical initial conditions 
of {x, p, T} and started off from the same initial ran- 
dom configuration. As an illustration, Figure [T] shows 
two corresponding structures obtained with and without 
Coulomb. 

The effect of the screened potential can also be ob- 
served on the fragment size multiplicity. Figure [2] shows 
the distribution of cluster sizes observed in 200 con- 
figurations of asymmetric matter (x = 0.3) at density 
p = 0.015 /m~ 3 and temperatures T = 0.1 MeV and 
T = 1.0 MeV. The figures corresponds to configurations 
with (cases A and C) and without (B and D) the Coulomb 
interaction. The increase of the number of fragments of 
sizes A > 500 as a goes from a — 1 to a = underline 
the role the Coulomb force has on determining the mass 
distribution. 



6 



The change of the inner structure can be quantified 
through the use of the radial distribution function, g(r). 
An example of this effect is shown in Figure |3] which 
shows this function calculated with the structures ob- 
tained with symmetric nuclear matter (x = 0.5) at den- 
sity p = 0.072 fm~ 3 and temperature T = 0.1 MeV and 
for varying strengths of the Coulomb potential, namely 
a = 1 (full Coulomb), 0.8, 0.2, and (without Coulomb). 
As it can be easily observed, the proton-electron gas in- 
teraction does not affect the nearest neighbor distances, 
although the magnitude of the function for second near- 
est neighbors increases as Coulomb weakens. 

The isotopic content on the cluster is also affected by 
the presence of the proton-electron gas interaction. Fig- 
ure H] shows the x content of the clusters formed in 200 
configurations of asymmetric matter (x = 0.3) at density 
p = 0.015 fm~ 3 and temperature T = 1.0 MeV . Again, 
the top panel corresponds to the case with full Coulomb 
and the bottom one to the case without such interaction. 
Clearly visible is the enhancement up to x w 0.4 that 
would occur without Coulomb; also noticeable is the in- 
crease of the maximum fragments sizes which grow up to 
A w 800 without the electric interaction. 

As suspected in [l^|, in this study we find that the 
interactions between protons and electrons are responsi- 
ble for the rearrangement of protons, although we find 
this to be of varying degrees. Figure [5] shows the RMS 
displacement of nucleons as a function of the evolution 
of the stationary ergodic process. Each curve is made of 
200 points, each of which represents the average displace- 
ment (with respect to the original configuration) of the 
2000 nucleons used in the x = 0.5 cases shown. For each 
case tracked (p = 0.015 and 0.72 /to" 3 at T = 0.1 and 
1.0 MeV), the displacements obtained with and without 
Coulomb are presented. Curiously, at low temperatures 
Coulomb seems to enhance the nucleon mobility much 
more than at high temperatures. 

A similar effect is observed in the persistence, which 
is the tendency of nucleons in a given fragment to re- 
main in the same cluster. Figure [6] shows the persistence 
with and without Coulomb as it evolves through 200 con- 
figurations of the stationary ergodic process, measured 
with respect to the original configuration of 3000 nucle- 
ons with x = 0.3, p = 0.015 and T = 1.0 MeV. As seen 
in the mobility, Coulomb enhances the transfer of nucle- 
ons thus decreasing the persistence as shown by the two 
curves. 

The effect of Coulomb on the shape of the structures 
formed can be quantified using the Euler Characteris- 
tic and Mean Curvature introduced before. Following 
the procedure explained in more detail elsewhere [l3j . 
digitized polyhedra were constructed for each of the 200 
nuclear structures obtained with and without Coulomb 
at each value of {x, p, T} considered. The corresponding 
values of the mean curvature and Euler numbers were 
calculated and, although the configurations correspond- 
ing to the same set of conditions were not identical, their 
curvature and Euler numbers showed robust average val- 
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FIG. 7: Average values of the Curvature and Euler num- 
bers of the structures listed in Table [Til circles correspond 
to structures with Coulomb and squares to structures with- 
out Coulomb, arrows indicate the average displacement of the 
structures as a goes from 1 to 0. 



ues. Table ITT1 shows these values and Figure [7] their loca- 
tion on the Curvature-Euler plane as well as their changes 
in position as the Coulomb strength is diminished. The 
standard deviations of the points are, roughly speaking, 
of the size of the points used in the plots. 



TABLE II: Classification Curvature - Euler 





P 


T 


a = 


1 


a = 







(fm~ 3 ) 


(MeV) 


Curvature Euler 


Curvature Euler 


A 


0.072 


1.0 


-1740 


89 


-24 


4.7 


B 


0.072 


0.1 


-382 


-15.8 


29.4 


-1.9 


C 


0.015 


0.1 


1771 


8.9 


514 


3.5 


D 


0.015 


1.0 


2212 


-1.7 


273 


8.6 


E 


0.015 


1.0 


9090 


-612 


11011 


-243 


F 


0.015 


0.1 


10804 


454 


11436 


159 


G 


0.072 


1.0 


-2969 


447 


-4512 


602 


H 


0.072 


0.1 


-4485 


798 


-5914 


933 
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A. Analysis 

The results obtained in this study allow us to under- 
stand the effect the proton-electron gas interaction has 
on the formation of the nuclear pasta. Perhaps the most 
interesting result of this investigation is the fact that the 
characteristic topological structures of the nuclear pasta 
-which are usually associated with the presence of the 
competing interactions of the long range Coulomb poten- 
tial and the short range nuclear force- exist even when 
there is no Coulomb (i.e. when a — > 0). We believe this 
to be the result of the competition between attractive 
and repulsive nuclear forces at very low temperatures; in 
our case this is due to the attraction of the V np and the 
repulsion of the V nn and V pp which operate at medium 
ranges; incidentally, we find the same effect also in other 
potentials [l9j as will be reported in a future communi- 
cation. 

Thus the role of the Coulomb force is to emphasize such 
competition between attractive-repulsive forces. The ex- 
amples of structures formed with and without Coulomb 
(Figure []} suggest that the electric forces tend to dis- 
tribute matter more and form less-compact objects. This 
behavior is in agreement with the reduction of the frag- 
ment size multiplicity undoubtedly induced by Coulomb 
as illustrated in Figures [5] and @] 

At a microscopic scale, Figure [3] shows that varying the 
strength of the Coulomb interaction does not change the 
inter-particle distance, but certainly decreases a bit the 
x content of the fragments as can be seen in Figure |U In 
agreement with this, the increment of nucleon mobility 
produced by Coulomb (cf. Figure [5]) gets reflected in a 
reduction of the persistence (cf. Figure [6]). 

With respect to the modification of the shapes due to 
the electric interaction, we find that Coulomb tends to 
reduce both the curvature and Euler number of the struc- 
tures. Indeed, Coulomb appears to "compactify" sym- 
metric matter into spaghetti and lasagna type structures 
as shown clearly in Figure [7] for symmetric matter (top 
panel); of course, the effect becomes blurrier in systems 
with smaller charge densities, such as the asymmetric 
matter case of x = 0.3 (bottom panel). 

IV. CONCLUDING REMARKS 

The effect of the proton-electron gas interaction in 
symmetric and neutron-rich matter was studied at den- 



sities and temperatures of interest to neutron star crusts 
by varying the interaction strength. Its effect on the 
fragment size multiplicity, the inter-particle distance, the 
isospin content of the clusters, the nucleon mobility, and 
on the modification of the topological shape due to the 
Coulomb interaction was studied. 

The most general result is the existence of the nu- 
clear pasta structures even without the presence of the 
Coulomb interaction due to the interplay between the 
attractive and repulsive parts of the nuclear forces. In 
general, the proton-electron Coulomb interaction tends 
to distribute matter more, form less-compact objects, in- 
crease nucleon mobility, reduce the persistence and the 
fragment size multiplicity, decrease the x content of the 
clusters, but does not change the inter-particlc distances. 
All of this gets reflected on a modification of the shapes 
of the nuclear structures formed by a reduction of their 
curvature and Euler numbers; effect that is more pro- 
nounced in symmetric than asymmetric matter. 

In conclusion, we find that the electron gas screen- 
ing is of upmost importance in the formation of the nu- 
clear pasta in neutron star crust environments. In a fu- 
ture work we will extend the present study to investigate 
the transition from a uniform nuclear phase to the pasta 
structures which is bound to occur at very low temper- 
atures, as well as its dependence on the isotopic content 
and compressibility of the nuclear matter using both ver- 
sions of the Padharipande potential as well as other clas- 
sical potentials [ljl |4(| ■ 
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